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Abstract 

We show the propagation of regularity, uniformly in time, for the scaled solutions of the 
inelastic Maxwell model for small inelasticity. This result together with the weak convergence 
towards the homogenous cooling state present in the literature implies the strong convergence 
in Sobolev norms and in the L 1 norm towards it depending on the regularity of the initial data. 
The strategy of the proof is based on a precise control of the growth of the Fisher information 
for the inelastic Boltzmann equation. Moreover, as an application we obtain a bound in the L 1 
distance between the homogeneous cooling state and the corresponding Maxwcllian distribution 
vanishing as the inelasticity goes to zero. 

Mathematics Subject Classification Numbers: 82C40, 35B40 



1 Introduction 

This paper concerns the regularity properties of solutions of the spatially homogeneous Boltzmann 
equation for Maxwellian molecules in R 3 with inelastic collisions, introduced in [5]. This equation 
describes the evolution of the distribution of the velocities in a collection of particles as they interact 
through inelastic binary collisions. Let f{v , t) be the probability density for the velocity of a particle 
chosen randomly from the collection at time t. Let (p be any bounded and continuous function on 
IR 3 . Then the equation under investigation is given, in weak form, by 

^(f,<P) = (Qe(f,f),<p) (LI) 
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partially supported by the DGI-MEC (Spain) FEDER-project MTM2005-08024. MCC is partially supported by FCT 
POCI/MAT/61931. 
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where (/,</?) denotes f R3 f(v,t)ip(v)dv, and where 

If (' f ~ ( v — w 



(Q*U,f)M 



47T 



f Bin- 

52 V \v — w\ 



f(v)f(w)[ip(v / ) — ip(v)]dndvdw (1-2) 



where n is a unit vector in S 2 , dn is the uniform measure on S 2 with total mass Air, with 1/ the 
post collisional velocity given by 

1 + e, 



-((t> — w) ■ n)n , 



(1.3) 



with < e < 1, and where B is a positive, integrable even function on [—1,1]. Because of the 
integrability of B, we can separate the collision operator in the gain and loss terms, Q e (f, f) '■= 

Qi(f, /)-«(/,/) with 

«(/,/) = (J s2 B (» • ^) dn) / = Q £ *) /. 

The function i? gives the rate at which the various kinematically possible collisions happen, and 
the tilde is present because later on we shall also consider another rate function B corresponding 
to another parameterization of the kinematically possible collisions. 

The parameter e is the restitution coefficient. For e < 1, the collisions are inelastic, and energy 
is dissipated in each collision. In this case, the collisions are not reversible. This is a crucial 
difference with the elastic theory in which there is a complete time reversal symmetry between the 



pre and post collisional velocities. It is for this reason that we have written (1.1 ) in weak form, and 



not because of any difficulty in constructing strong solutions: It is just that to write the equation 
down in strong form, we need a parameterization of the possible pairs of precollisional velocities 
(v*, w*) that can result in the pair of post collisional velocities (y, w). We have detailed the strong 
formulation in the different parameterizations in the Appendix, see O [18] for particular cases, but 
for the purposes of the introduction, the weak form specifies the equation under consideration well 
enough for us to proceed with the description of the particular issues with which we are concerned, 
and the results we obtain. 



The reason we require B to be even is that the post collisional velocity v' defined in ( 1.3 ) depends 
on n quadratically, and thus is unchanged under the substitution n — > — n. For this reason, we 
may freely assume B to be even, and we do so in what follows. 

The first thing to notice about the equation is that the first moment of / is conserved. Indeed, 
I 3 , let if(v) = a ■ v. Then we have from (1.1 ) and (1.3) that 



for any a £ 



d 



1 + e 
1 + e 



p , v-w 
B n 
S 2 \ \v — w 



s 2 B(s)ds 



f(v)f(w)[((v — w) ■ n)(a ■ n)]dvdwdn 
f(v)f(w)a ■ (v — w)dvdw = . 



Indeed, as detailed in the appendix, the companion formula to (1.3), giving the other post collisional 
velocity w', is 

w' = w H — — ((t> — w) ■ n)n . 
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Thus, in each individual collision (v,w) — ► (v',w'), the total momentum v + w is conserved, and 
this certainly ensures that the first moment of / is conserved. In any case, on account of the 
computation just made, we may as well assume that our initial data fo satisfies 



vf (v)dv = 



Then of course we shall have 



vf(v,t)dv = (1.4) 

for all t > 0. While momentum is conserved, energy is dissipated, as we have indicated above. We 
now calculate the rate of this dissipation: Take tp(v) = \v\ 2 and then note that from (1.3) 

ip(v ) = r 



(1 + e) 2 

' 2 (1 + e)((v - w) ■ n)(v ■ n) + ^ — ^-((v ~ w ) " n ) 2 • 



In this case, using the abreviated notation u = v — w, we have 

■(l + ef 



1 

47T 
1 

2 



S 2 V \U\ 



s 2 B(s)ds 



-i 



s 2 B(s)ds 



Hv)f(w) 



■<J,<P) ■ 
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(1 + e) 2 



(u ■ n) 2 — (1 + e){u ■ n)(v ■ n) 
\u\ 2 — (1 + e)u ■ v I dvdw 



dvdwdn 



That is, with the positive constant E defined by 



E 



s 2 B(s)ds 



we have every solution of (1.1) with initial data fo satisfying ([TJ satisfies 



\v\ 2 f(v, t)dv 



-2E 



\v\ 2 f(v , t)dv . 



(1.5) 



(1.6) 



This implies that f(v, t)dv tends to a point mass at v = as t tends to infinity. It is natural to 
enquire into precise nature of this collapse to a point mass. In previous works [6j |9j El [71 [H [10], it 
has been shown that if one rescales f(v,t) to keep the variance (i.e., temperature) constant, then 
the rescaled density tends to a particular equilibrium state, known as the homogeneous cooling 
state. That is, if we define the probability density g(v, t) by 



g(v,t) 



-3Et 



-Et„ 



(1.7) 



\v\ 2 g(v, t)dv 



\v \ 2 g(v , 0)dv 



\v\ 2 f (v)dv 



for all t, and there is a density g^ such that 

lira g(v,t) = g^v) . 



t— too 



(1.8) 



The convergence in (1.8), part of the so-called Ernst-Brito conjecture [HI (TTJ El [9] , has so far been 
shown in certain weak norms that we shall introduce shortly, see [13] for a review and [91 [2] for 
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the proofs. Our goal in this paper is to prove that g(v, t) is regular in v, uniformly in t. This is 
reasonable to expect since it is was proved by Bobylev and Cercignani that goo(v) itself is quite 
regular [61 Theorem 7.1]. However, it is clear from the fact that f(v,t)dv tends to a point mass that 
the norm of /(•, t) must diverge with t in every norm that would imply smoothness of /(•, t). While 
the rescaling may well lower such norms, one needs very precise estimates on the loss of regularity 
to avoid having them overwhelm whatever one gains from the rescaling. Notice in particular that 
the rescaling does nothing to improve regularity. 

To investigate the long time behavior of g(v,t), we write down its evolution equation, which 



of course is obtained from (1.1) through (1.7). In working out the equation, we make use of the 



dilation invariance of the collision integral Q e (f,f): For any density /, test function ip and any 
A > 0, define 



f W (v) = X 3 f(Xv) 



and 



Then one easily sees from (1.2) that 



{Qe(f W ,f W ) t V>) = {Qe(f,f),<P W ) 



(1.9) 



;i.io) 



Then, for any test function </?, 



d {/ (-PC-**))) 
dt 



-Ee- Et (f,v • (Vc^ cxp (-^) + (Q(/,/)^ (cxp( -^») 
-E(g,(vVp)) + (Q(f,f),^- Et ») 
-E(g,(vVtp)) + (Q(g,g),ip) 



where we have used (1.10) in the last line. Thus, our evolution equation for g is, in weak form, 



dt 



(g,<p) = -E(g,{vVcp)) + {Q{g,g),ip) . 



a. in 



There are other ways, physically and mathematically different, to control the tempera- 
ture/variance: If the particles are in contact with an appropriate heat bath, this will add a thermal 
regularization to the evolution equation for /. This thermal bath can be modelled by stochastic 
heating, i.e., 

^ = Qe(f,f) + A v f 

or by a thermalized bath of particles, adding a linear Boltzmann type operator. In these two cases, 
global regularity estimates for solutions have been obtained, see [U [23] . However, as the first 



order anti-drift term in (1.11) does not induce a priori any regularization, the problem of proving 



global regularity estimate for solutions of ( |1.11[ ) is more challenging. 

The Fisher information plays a crucial role in our investigation of regularity. For any probability 
density / on M 3 , the Fisher information, 1(f), is defined by 



1(f) 



V v r f(vj\ 2 dv= \Vlnf(v)\ 2 f(v)dv 



whenever the distributional gradient of v7 is square integrable, and it is defined to be infinite 
otherwise. It has been shown by Villani that in case e = 1; i.e, for elastic collisions, the Fisher 
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information is non increasing in time. This is a basic propagation of regularity result that is the 
starting point of our investigation of the inelastic case. 

For solutions of ( |1.11[ ), the Fisher information will not be bounded uniformly in time. Indeed, 
the Fisher information has simple scaling properties: If is defined in terms of / and A as in 
(1.9), one easily computes 

J(/W) = X 2 I(f) . (1.12) 



Therefore, with g(v,t) defined in terms of f(v,t) through the scaling relation (1.7), we have 

I(g(;t))=e~ 2Et I(f(;t)) . (1.13) 



The exponentially decreasing factor e is good, but notice from (1.5) that it depends only on B 
and the restitution coefficient e, and not on the initial data. For some initial data, /(/(•, t)) will 
grow faster than this rate, and thus I(g(-,t)) will grow exponentially. Nonetheless, we shall be able 
to prove that its growth is not too bad, at least for e not too far from 1. 

1.1 THEOREM. For any solution g(v,t) of ( |1.11[ ) ; we have a bound on the Fisher information 

I(g(;t)) < e ^-^ + e + l^)/(8e^)-2E]t I{g{ ^ 0)) 



where E is the constant defined in (1.5). While the exponent is always positive, no matter how B 
is chosen, it does always vanish in the limit e— ► 1. 

Theorem |1.1| is proved in Section 2. Our main goal in the next part of the paper is to obtain a 
tiny uniform-in-time propagation of regularity result of the type: 

1.2 THEOREM. For any < 5 < 1, there is a computable positive constant C, such that for 
any solution g of ( |1.11[ ) corresponding to the initial value fo with unit mass, zero mean velocity, 
\v\ 2+a fo G L 1 ^ 3 ) with < a < 1 and I(f ) < do, then 

\\\v\ S 9(v)h°°(R3)<C, (1.14) 

for all t > 0, being e close enough to 1. 

This strategy precisely coincides with the open problem left in |13| for strong convergence to 
homogeneous cooling states and applied in the case of the thermalized bath of particles, adding a 
linear Boltzmann type operator, see [HI Subsubsection 7.2.4]. 

To prove the convergence in strong norms towards the homogeneous cooling state, we will need 
more; we will need the propagation of regularity in Sobolev spaces of high degree: 

\\9\\lr { ^ = I ' M*\m?*n (i-iB) 



with r > 0. However, there is a well-developed machinery |12| [T3] for showing that whenever the 
equation propagates a tiny degree of regularity, this implies the equation propagates regularity of 
any degree. Therefore, the main problem to be solved is to prove (1.14), uniformly in time for 
which there are no standard arguments. 

Then, using the regularity in high Sobolev spaces, we can parley the weak convergence in (1.8) 



into convergence in all Sobolev norms, Theorem 3.9, and strong L 1 convergence at an explicit 
exponential rate for a certain class of initial data. This is the objective of Section 3 and the main 
result is summarized as: 
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1.3 THEOREM. Given the solution g of corresponding to the initial probability distribution 

function /o E i? r (IR 3 ), with r > 0, of zero mean velocity such that |f| 4 /o E L 1 (M 3 ) and I(fo) < oo. 
Then, for e close to 1, the solution g(t, v) of ( |1.11[ ) converges strongly in L 1 with an exponential 
rate towards the homogenous cooling state, i.e., there exist positive constants C and 7' explicitly 
computable such that 

\\g(f) - 9oo||li(r3) < Ce" 7 '* 

for all t > 0. 

Finally, we can study the small inelasticity limit of the sequence of homogeneous cooling states 
showing an L 1 convergence towards the Maxwellian distribution with zero mean velocity and tem- 
perature fixed by the initial data as e — » 1 with an explicit speed in terms of the inelasticity 
parameter. Section 4 is devoted to this small inelasticity limit in strong norms. Finally, as an- 
nounced above, the appendix is aimed at a detailed description of the relations between the different 
parameterizations of the collision mechanism that we have written for non necessarily Maxwellian 
type collision kernels. 



2 Fisher Information bounds 

Villani [23] has proved that for Maxwellian molecules and elastic collisions, the Fisher information 
does not increase. A special case of this, namely with B constant, had been treated earlier by 
Carlen and Carvalho using the reflection parameterization [11]. Villani's analysis is based on the 
a representation, and has the advantage it allows an arbitrary rate function B. All of these results 
use the strong formulation of the collision operator, the passage from the weak form to the strong 
form is merely a complicated change of variables that we detail in the appendix. The main formulas 
we will use in this section are related to the strong formulation in the cr-representation, 



4tt 



with u = v — w, k = u/\u\ 



r, ( , P , (l + e 2 )*-(l-e 2 ) \ V2 1 B (V (1 - -)/2 



(1 + e 2) _ (1 _ e 2 )s J y (1 + g2) _ (1 _ e2)s e ■ 2y/(l-s)/2 

and the precollisional velocities are given by 

v + w 1 — e , . 1 + e. . 

v = — (v — w H : — \v — w\a 

2 4e v ; 4e 1 1 

v + w 1 — e, . 1 + e. . 

w = (v — w) \v — w \a 

2 4e v ' 4e 1 1 

The reader can understand now why we have avoided the strong formulation as long as we could. 



We emphasize that this operator coincides with the one defined in weak form below (1.2). Full 
details of the passage from one representation to the other are given in the appendix. 

We now start to adapt Villani's analysis to the inelastic case, and derive bounds on the growth 
of the Fisher information in terms of the restitution coefficient e. The crucial feature of these 
bounds on the growth is that they vanish as e tends towards 1. The main result of this section is: 
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2.1 THEOREM. For all probability densities f on M 3 , 



l(Q+(f,f))< 



1 + (1 - e) 



2 + e + 15e 



2\ 1 



'(/) 



TOi/i i/ie consequence that if f(v,t) is a solution of (1.1), we have 



(2.1) 



(2.2) 



/(/(•, t)) < e [( 1 - e )( 2+e+15e2 )/( 8e3 )]*i(/(- ) o)) . 

As an immediate consequence of this, we obtain the proof of Theorem 1 1.1} 
Proof of Theorem |1.1} Consider any rescaled solution g(v,t); i.e., any solution of ( |1.11[ ). By 



(1.13), any solution g(v,t) of ( 1.11| ) satisfies 

I(g(;t)) < e [(l-)(2+e+15e^)/( 8e 3)-2B] tj(5( . )0)) 



(2.3) 



where E is given by (1.5). Notice that while E depends on the particular choice of B, for any 
choice we have 

^ / , 1 + e 1 — e 
£<(l-e)— < — ■ 



Therefore, for any B, the exponent in (2.3) is at least 



( /2 + e + 15e 2 
(1_e)( 8^ X ) >0 ' 



for all < e < 1. 



□ 



While the exponent in Theorem 1.1 is always positive, it does vanish in the elastic limit, and 
that is what we shall need in the next sections. We begin by recalling several results: 

2.2 LEMMA. |23} Lemma 1] Let P^ denote the orthogonal projection onto the span of k, and let 
Pjt~ denote its orthogonal complement. Then for any difjerentiable rate function B, 

1 



V„ [B(k ■ a)} = ri B'(k ■ o-)P£<t 



where B' is the derivative of B. 



The proof of this Lemma is an elementary computation. Applying it with B = Bj, and defining 
F(v,w,a) = f{v*)f(w*) we have 



V„Q+(/,/) 



1 

4-7T 



1 



(Bj)'(k ■ a)P^aF{v, w, a) + B+{k ■ a)V v F(v, w, a) 



dwda . (2.4) 



The proof of the next Lemma is not so elementary, as done in [23], it is an ingenious integration 
by parts on the sphere. Later on, we will give a different proof of the formula resulting from this 
Lemma in a more direct way. 
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2.3 LEMMA. [23, Lemma 2] Using the notation of the previous lemma and also defining the 
linear transformation M a k on ^ 3 by 

M ai k{x) = (k • a)x — (k ■ x)a , 

we have that for any smooth function F(a) on S 2 , 

[ B'(k ■ a)P^aF{a)da = [ B{a ■ k)M a k [V a F(a)] da . 
Js 2 Js 2 



Using this Lemma on (2.4) we obtain 
1 



47T 



S 2 



1 



u 



B+(a ■ k)M aM [V a F(v, w, a)] + B+(k ■ a)V v F(v, w, a) 



dwda . (2.5) 



Next, using (A. 36) to evaluate the Jacobians, 



and 



V - = (^J (V«/)(«*) = ^-(V„/)(W) + • V v f)(v*)k 

V„ [f(w*)) = f 9 ^) (V„/)K) = ^ (V„/)K) - ^(a ■ V v f)(w*)k . 



Also from (A. 36) 



dv 



V a F(v,w,a) 



4e 



4e 



l + e 
4e 



u\[(V v f)(v*)f(w*)-f(v*)(V v f)(w*)} ■ 



Therefore, if we define the linear transformation P a ^ on M 3 by P a ,k{x) = (a • x)k + M a ^(x), we can 



rewrite (2.5) as 



where 



V„Q+(/, f) = ]~ f [ B+(a ■ k)G(v, w, a)dadw 

47r JS 12 7R3 



G(v,w,a) = /(«,*) ( + ) [(V„/)(0] 

+ /f<- J ) ( ^ - ) [(v,/)K)] 



(2.6) 



4e 



4e 



(2.7) 



Before proceeding further, let us give a simple, direct proof of formula (2.7) making use of the 
Fourier transform instead of Lemma 12.31 



Proof of (2.6)-(2.7): We start by recalling the formula of the Fourier representation of Qf(f, f) 
obtained in [5] , see previous works O H] and |15| [TB] for a review. It holds 

1 



with f] = 7]/\rj\ and 



Qt(f, f)(v) = r / B(fj- a)f(t, r,-)f(t, r)+) da 



n- = -^-(v -\vW), 



4 

3 - e 



V+ 



j] + 



l + e 



(2.8 



(2.9) 



\a = T) - rj. 
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Now, let us point out the following identity, left for the reader to check, 



vt \ 3e ~ 1 , 1 + e D / ^, 1 + e 



1 + e 1 - e 2 



4e 4e 4e 4e 4e 

Now, multiplying both sides by j=B(fj ■ a)f(t,r]_)f(t,ri + ) and integrating over the sphere, we get 

±- f B^-a)f(t,^)f(t,7 1+ )Z( V+ ,7 1 ^da = riQt{rf)(v) = [^v^,f)}(v) (2-10) 

since the integral of the last term is zero. In fact, since we are free to choose our coordinate system 
in the sphere, it is easy to see that 



S 2 



(ri ■a)r)J(rj,a)da = / \rj\a J{r), a 
JS 2 



da 



for any function J(rj,a). The desired formula (2.6)-(2.7) is just the inverse Fourier transform 
formula corresponding to (2.10). □ 



Now, starting from ( |2.7[ ), we can define H(v,w,cr) by 
H(v, w, a) = 2y/f(w*) (^fl + X -^P aM 



4e 4e 
1 + e 1 + e 



4e 



4e 



-Prr 



(v«v7)(«*) 
(v„>/7)(«0 



3e - 1 
2e 



yM(v v y/)(0 + ^//VT(v^)K) 



2e 



1 + e 
2e 



: Hx(v, w, a) + P a , k H 2 (v, w, a) 



(2.11) 



where the last line defines Hi(v,w,a) and H2(v,w,a), and G(v,w,a) = \J f(v*)f(w*)H(v, w, a) 
Thus, we get 



V,Q+(/,/) 



4?r 



5 , 



S+(cr • k)\/f(v*)f(w*)H(v, w, a)dadw . 



(2.12) 



1/2 



Therefore, by the Schwarz inequality 

\V v Qt(f,f)(v)\<(Qt(f,f)(v)) 1/2 (^[ [ B+(a-k)H 2 (v,w,a)dadw) (2.13) 

From here we obtain a bound of I(Qj(f, /)): Squaring both sides, and integrating in v we obtain 

HQtU,f))<T- I I I B+(a-k)\H{v,w,a)\ 2 dadwdv . (2.14) 
4vr J R3 J R 3 J S 2 

It remains to estimate the integral on the right in terms of /(/)■ This consists of the sum of three 
terms: 



Air 



S, 



B+(a ■ k^H^v^w^a^dadwdv 
B^io ■ k)\PcrfiH2(v, w, a)\ 2 dadwdv 



Sa= — 



47T y R s y R3 y 5 2 

7-111 B+(a-k)(H 1 -P Uik H 2 )(v,w,a)dadwdv . 

47r Jr3 y K 3 y 5 2 



(2-15) 
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Summarizing the discussion so far, we have: For any probability density / on M 3 , we have 

I(Qt(fJ))<S 1 + S 2 + 2S 3 , 



where the quantities on the right hand side are specified in (2.15 ). Our next lemma simplifies these 
expressions by a change of variables: 

2.4 LEMMA. For any probability density f on M 3 , we have 



Si 
S 2 



1 

47T 
1 

4tt 



S 2 



g2 



B(a ■ k^F^v.w^dadwdv 
B(a ■ k)\P a , jk ,F 2 (v,w)\ 2 do-dwdv 



(2.16) 



S 3 = ~ [ [ I B(a- to) ■ P a r, k ,F 2 )(v, w)dadwdv 

47r JR3 JR* J S 2 



where 



k' 



(1 - e)k + (1 + e)a 
■v/2(l + e 2 ) + 2(1 - e 2 )k ■ a 
'3e - 1 



a 



(1 + e)k + (1 - e)a 
v / 2(l + e 2 ) + 2(l-e 2 )A:-a 



and 



Fi(v,w) 



F 2 (v,w) 



1 + e 



. v v7(ii)(vv7)(«) + -^-v7W(vv7)H 



l + e 
2e 



V?H(vv7)(«) - vTMCVv^/)^: 



(2.17) 
(2.18) 

(2.19) 



Proof: In the expressions in (2.15), we are integrating over post collisional variables. We use the 



change of variables Theorem A.l which concerns the transformation C s , e {v, w, a) i— ► (v* ,w* ,a*) 



from post to pre collisional variables under the "swapping map"; i.e., for the sigma representation. 



Consulting Theorem A.l and the definition of H\ and H 2 in (2.11 ), we see that each of the integrands 



above can be written out in the longer form appearing in Theorem A.l e.g., 

\Hi(v,w,a)\ 2 = K[(v*,w*,a*),(v,w,a)} = K[C~}{v, w, a), (v, w, a)} 



Theorem A.l allows us to write this as an integral over 

K[(v, w, a),C Ste (v, w, a)} = K[(v, to, a), (v , w', a')} 



Doing this for each of the three integrals in (2.15), we obtain the stated formulas. 



□ 



Define the matrix A a ^ by 

■4<r,fe = Pa',k' ~ Pk,a ■ 

We shall now prove: 

2.5 LEMMA. For all a, k and z in M 3 such that \a\ = \k\ = I, 

\A a , k {z)\ < 2 ] — 6 \z\ . 

e 

In proving this Lemma, as well as for estimating S 2 , we shall make use of the following lemma 
of Villani: 
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2.6 LEMMA. |23j Lemma 4] For all a, k and z in R 3 such that \a\ 

\Pa,k( z )\ < N 

with equality if and only if a, k and z belong to the same plane. 



Proof of Lemma 2.5:, considering the formulas for k' and a' given in (A. 25) and (A. 26) respec 



tively, notice that as e — * 1, we have k! — > a and a' — > k, as we should, since in the elastic case, 



this is what the swapping map does. Therefore, using (A. 25) and (A. 26) we compute 

(1 - e)k + (e + 1 - y/2(l + e 2 ) + 2(1 - e 2 )k ■ a) a 



k'-a 



y/2(l + e 2 ) + 2(1 - e 2 )k ■ a 



and 



z , (1 - e)a + (e + 1 - ygg + e 2 ) + 2(1 - e 2 )k • a)k 

^2(1 + e 2 ) + 2(1 -e 2 )k-a 

Using the elementary estimates 

2e < y^2(l + e 2 ) + 2(1 — e 2 )k ■ a < 2 , 

we easily find that 



\k' -a\ < 



1 - e 



and 



W -k\ < 



1 - e 



e e 
Now, notice that P a > ^ = + P a >-k,k' + Pk,k'-a- This means that A U)k = P a >- k ,k> + Pk,k'-a, and 
now the result follows from Lemma 2.6 and the triangle inequality. □ 



Now we are ready to estimate Si, S2, and S3 in terms of /(/), and prove the main result of this 
section. 



Proof of Theorem 2.1 : First of all, notice that by Lemma 2.6 

S2 < S2 

where 

S 2 = — / / B((T- k)\F 2 (v,w)\ 2 dadwdv 
4vr J K 3 J K 3 J S 2 



(2.20) 
(2.21) 



Next, we have 

Sz = ^- [ [ [ B(a-k)(Fi-P ki(T F 2 ){v,w,a)dadwdv 

^ JR3 J R 3 J S 2 

+ 7- f f [ B{a-k){F 1 -A ak F 2 ){v,w,a)dadwdv. (2.22) 

4tt J K 3 J K 3 J S 2 



Since B is even in a, P k ,a is odd in a, the first integral is zero. Then, by Lemma 2.5 and the 
Schwarz inequality, 

S3 < 2 l —Js~J 2 < — (Si + S 2 ) • (2.23) 
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Notice that this vanishes in the elastic limit; this is the key point discovered by Villani in the elastic 
case. Finally, we need to estimate Si and S 2 - For this, notice that 



\Fi(v,w)\' 



+ 2 



3e 2 + 2e - 1 



4e 2 



l + e 
2e 



and 

|F 2 (w,u;)| 2 = 
Thus, 

\Fi(v,w)\ 2 + \F 2 (v,w)\ 2 



/HI V vTMl 2 + /(«)|V vTHI 2 - 2y/f(v)f(w)Vy/f(v) ■ Vi/Jivj 



5e 2 - 2e + 1 



+ 



2e 2 
e 2 -l 



/H|v v / 7Ml 2 + 



e 2 + 2e + 1 
2^ 



/(^IvTThP 



y/f{v)f{w)Vy/J{v) ■ V vTM • 



Now using the fact that, with our chosen normalization, /« 2 • cr)da = 4ir, together with the 
Schwarz inequality and the definition of 1(f), we have 

7e 2 + 1 



8e 2 



Combining this with (2.23), we obtain, 



Si + S 2 + 25 3 < 1 + 2 



1 - e\ 7e 2 + 1 



/(/) = 1 + (1 - e) 



2 + e + 15e 2 



1(f) ■ 



This proves Theorem 2.1 



□ 



3 Propagation of regularity 



The next lemma relates the Fisher information bound to an L°°-Fourier bound, similar arguments 
were used in |21| [13]. Nevertheless, we include its idea for completeness. 

3.1 LEMMA. For any probability density g on M 3 , there is a constant C such that 

\\\v\9(v)\\loo(r 3) <CI(g) 1 ' 2 . 

Proof: Let h = yjg. Then, the Fourier transform of g can be written as the convolution of h 
with itself, g(ri) = (h * h)(rj). Now, the boundedness of the Fisher information of g implies that 
h= ^ge H 1 ^ 3 ), and thus 



\v\\d(r])\ = \v\ 



h(rj — r]*)h(r]#) dr]* 



< 



(\rj - 77*| + M) \h(rj - rj*)\\h(ri*)\dri* 



<2 



\ V \ 2 \h(v)\ 2 drj 



1/2 



\h( v )\ 2 d V 



1/2 
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giving the desired result. 



□ 



From now on, we will restrict our attention to the most relevant case in the literature in which 
B is constant and due to normalization B = 1 or equivalently B(s) = 2\s\, see appendix and ( |A.19 ). 
In this case, we remind from (1.5) that E = (1 — e 2 )/8. We shall combine Theorem 1 1 . 1 1 with the 
following result, due to Bobylev, Cercignani and Toscani [9j and Bisi, Carrillo and Toscani [2] in 
this form, see also previous results [6], which gives the uniform weak norm control. We need some 
notation, for any < a < 1, let us consider 



A + a 

and A 2 (a,e) = 1 - Ai(a,e) - E (2 + a). 



1 + e 



2+a 



l__ e N4+a 



+ 



i-(¥) 2 



3.2 THEOREM. [2J Theorem 4.6, Remark 4.7] For any solution g(v, t) of ( |1.11[ ) with constant 
collision frequency B = 1, corresponding to the initial value fo with unit mass, zero mean velocity 
such that \v\ 2+a fo £ L 1 (IR 3 ) with < a < 1, there exist positive constants C(fo) and 7(0, e) such 
that 

d 2 (g(t),g 00 



sup l?(7? ' t) ,7 oo(r?)l < Ce* 



for all t > with 



mm 



2 (3-e)(l + e) 
M{oL,e), 



2 + a 



Let us remark that 



j(a, e) — > 7* := min 



2a 



(2 + a)(4 + a)' 2 



2 

15 



(3.1) 



as e — > 1 and a — ► 1 respectively. Combining Theorems 1 1 . 1 1 and |3 . 2 1 and Lemma 3.1 we shall prove 
one of our main results, Theorem |1.2[ whose statement we now make more precise. 



3.3 THEOREM. For any < 5 < 1, there are computable positive constants C, 7', such that for 
any solution g of corresponding to the initial value fo with unit mass, zero mean velocity, 

\v\ 2+a f € L 1 ^ 3 ) with < a < 1 and I(f Q ) < 00, then 

sup <Ce-T'*+ sup \v\ 5 \doc(v,t)\ =Ce-^ i + C 00 



/or all t > 0, £>em<7 e dose enough to 1. 



Proof: Pick some R > 0. By Lemma 3.1, for all r/ with |r/| > R, and all 5 < 1, 



where we used Theorem 1 1 . 1 1 and 

;i-e)(2 + e + 15e 2 ) 



)i 



ci(e) 



16e 3 



E . 
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On the other hand, for \r]\ < R, we have 

< ^ l^)-f^l +H ^l ?0o(77 )l 

< R^Ce-^ 1 + \r ] \ s \g 00 {r ] )\ . 
Combining estimates, we have that for all rj, 

\v\ S \d(il)\ < R 5 - l Ce tc ^ + R 5+2 Ce-^ 1 + M 5 |?oo(r/)| • 

We now minimize in R. Up to a constant multiple, the optimal choice is R = e^ 1 ^^ 7 ^' 6 ^/ 3 . This 
results m 

M*r/M^ ^ Cl (e)(J + 2) + ((5-l)7(a,e) ^ g 

15(77)1 < Cexp I t I + |r?| |</oo(»/)| • 

Choosing 5 < 1 we see that for e sufficiently close to 1, so that ci(e) is sufficiently close to and 



j(a,e) ~ 7* > 0, see (3.1), the exponent is negative. Finally, taking into account the regularity 
obtained by Bobylev and Cercignani for the homogeneous cooling state c/oo in Theorem 5.3], we 
deduce 

exp{-|r7| 2 } < |Soo(M)| < exp{-M}(l + \ V \) , 
from which Coo < oo. □ 

Now, let us proceed to write the evolution of Sobolev-type norms for our model. Since moments 
in Fourier space will have simpler relations, we shall use the homogeneous Sobolev quantities, with 
r > 0, defined in (1.15). Its evolution for solutions of ( |1.11[ ) is given by 



d_ 
dt 



\vr\9iv)rdv=2 I _J \vr9(v-)9(v+W(ri)dtrdTi 



2 ( |r / | 2 n?(r ? )| 2 (ir ? -2(2r + 3) / \ V \ 2r \g(v)\ 2 dr), (3.2) 



where z c is the complex conjugate of z. Let us start by estimating the contribution of the first 
term. We need to estimate the regularity contribution of Qf(g, g) and for this, we make use of the 
estimate of 

II M*iK*> *7) l|2,°°(R3) 



obtained in Lemma 3.1 In fact the situation is quite similar to [13, Subsubsection 7.2.4] and (T3| 
Lemma 7.13] in the case of thermalization by a bath of particles, adding a linear Boltzmann type 
operator. We will make use of the following lemma of Carrillo and Toscani. 

3.4 LEMMA. [131 Lemma 7.13, Proposition 7.30] Let g G H r {M?) and a probability density, then 
if 

III??I 5 9(^)IIl->(k3) < 00 
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holds with < 5 < 1 and r > f, i/ien 



|2r~ 



5-' 



g(v~)g(v+)g c (v) do- &r\ 



< C(r, e) \\\v\ S g(v)\\L°°(R3) WgW'jjr-s/: 



Here, the constant C degenerates as e — > 1 as C(r, e) ~ \ 2 4 



r _ 1 « 
4 



Taking into account Lemmas 3.1 Theorem|3.3|and|3.4[ we deduce from the evolution of Sobolev- 



type norms in (3.2) that 



— II II 2 
dt m ^ 



< DxWgf^^ -A(r + 2)\\g\\% r 



(3.3) 



with 



Di := C(r,e) sup |||ry| g(t, < oo. 

t>o 



We finally use standard Nash-type inequalities, see for instance |13[ Lemma 7.14]. 

L , . u _ 1 jrfjn>3\ i „ j,„z,„-7„-+„, A „•+ a-h, „ ^ S 

H r - S / 2 (R 3 ) and 



3.5 LEMMA. Let g G H r (M. 3 ) and a probability density with r > \, < 6 < 1, then g G 



with 



f \ (2r+3)/(2r+3-<5) 

\\g\\m(R3) - C r,S ^l|3llHr-«/2(R3) J (3-4) 

/ x \ 2/(2r+3-<5) / 2r + 3 _ 6 ^ (2r+3)/(2r+3-<5) 

^"U^J ( 2r + 3 



The previous lemma allows us to obtain the inequality 
d 



dt 



\g\\% 



< Do 



I II 2 
\9\\ht 



with D2 easily obtained from above and 6 = 
achieve one of the main theorems of our work. 



-4(r + 2)|| 5 ||| r(R3) (3.5) 
(2r + 3 — 5)/(2r + 3) < 1. As a consequence, we 



3.6 THEOREM. Given the solution g of (1.11) corresponding to the initial value fo G i/ r (IR 3 ), 
with r > 0, of unit mass, zero mean velocity such that \v\ 2+a fo £ L 1 (M 3 ) with < a < 1 and 
I(fo) < 00 • Then, for e close to 1, the solution g(t,v) of ( |1.11[ ) is bounded in H r (K 3 ), and there is 
a universal constant A so that, for all t > 0, 



IIs(*)IIh--(r3) < maxjll/oll^r-^),^} • 



In particular, the stationary solution or homogenous cooling profile g^ to (1.11 ) belongs to -ff°°(IR 3 ). 



3.7 Remark. Let /o G if r (M 3 ), with r > of unit mass, zero mean velocity such that |f| 2+a /o G 
L 1 (M 3 ) with < a < 1 and I(fo) < 00, 6e any initial datum for equation ( |1.11[ ). Previous theorem 



together with the Nash inequality in Lemma 3.5 implies that, for e close to 1, the solution g(t, v) 
of (1.11) is bounded in L 2 (M 3 ), and there is a universal constant C2 so that, for all t > 0, 



L 2 (R3) < max || f 



i3/(3+2r) 



c 2 }. 
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Let us point out that the previous propagation of smoothness results are true for any value of 
e for which a uniform in time estimate of || \r]\ s g(t, rf) [[^oo^) is available, which in our case is given 
by the values of e for which the estimate in Theorem |1.1| is satisfied with 7' > 0. 

Finally, using the strategy already introduced in [12] and used in inelastic models in [1] [13] , see 
also [24 , we can obtain the convergence in L . The first ingredient is an interpolation inequality 
that allows to control distances in arbitrary Sobolev norms and in 1? by using the propagation of 
smoothness and the convergence result in Theorem 3.2 in [2, Theorem 4.6, Remark 4.7]. 

3.8 PROPOSITION. [12, Theorem 4.1] Let s > 0, and fix > 0, < (5 2 < 1 be given. Then 



11/ - 9\\Hs m < C((3 U (3 2 ) d 2 (f, g)W> min (||/ - <?||^ 1(R3) , ||/ 



9\\H r 



with 



and 



s + 2(l-/? 2 ) 
C{f3 u f3 2 )-- 



r 2 



2s + (7 + &)(!- ft) 



4ir 



(1 + 3/Pi 



2/? 2 

l-/3 2 



sup 



\m-g(v)\ 



The previous result implies immediately convergence in strong norms: 



3.9 THEOREM. Let g be the solution of (|1.11|) corresponding to the initial probability distribution 
function /o E H r+e ' 



), with r > and e > 0, of zero mean velocity such that /(/o) < 00. Then, for 
e close to 1, the solution g(t, v) of ( |1.11[ ) converges strongly in H r with an exponential rate towards 
the homogenous cooling state, i.e., there exist positive constants C and 7 explicitly computable such 
that 

\\g(t)- goo \\ Arf ^<Ce-v 



for all t>0. 



Let us point out that the exponential rate 7 can be computed as 7 = (1 — (3 2 )"f for any choice 
of Pi > 0, < (3 2 < 1 such that max(ri, r 2 ) < r + e. Moreover, the uniform control of moments for 
the probability measure yields control of the distance in L 1 . 

3.10 LEMMA. [T2 Theorem 4.2] Let f G L 1 n L 2 (M 3 ) with \v\ 2p f G L 1 (IR 3 ), then, for all p>0, 



\f(v)\dv<C(p) 



\f(v)\ 2 dv 



\ 2p/(3+4p) 



with 



C{p) 



3 \ 4p/(3+4p) 

Ap) 



4p\ 

3; 



3/(3+4p)' 



\ 3/(3+4p) 

\v\ 2p \f(v)\dv) 



/ 47r \ 2p/(3+4p) 
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Let us recall what is known about tails of the homogeneous cooling state. One very interesting 
property is that not all moments of g^ are bounded and the threshold moment depends on the 
restitution coefficient e. This was proved by [6], see [8] for generalizations. In particular, the 
fourth moment of g^ is bounded for all restitution coefficients. 

The next ingredient we need to pass from 1? to L 1 convergence is the uniform-in-time propa- 
gation of the fourth moment for solutions for any value of the restitution coefficient as obtained in 



TUl Appendix] from which the main Theorem 1.3 immediately follows. 



The exact value of the constant 7' depends on the value of r in the hypotheses of Theorem 1.3 
In fact, by taking s = in Proposition 3.8 we get that the decay in 1? will be given by a constant 



7 = 7(1 — P2) for any choice of 0\ > 0, < 02 < 1 such that 

^ 2(1 -a) (7+ftKi-ftn 

Once, we have this decay of the I? norm, the previous lemma finally gives the value 7' = ^-7 



4 Small inelasticity limit of HCS 

As a further application of the results proven in the previous sections, we study the small inelasticity 
limit of the homogeneous cooling states and prove, as one might expect, that as e — > 1 then 
the homogeneous cooling state converges towards the corresponding Maxwellian in strong norms. 
Previous results in this direction were done in the asymptotic expansion in Fourier for the self- 
similar solution, see p2 Subsection 6.1]. 

Let us fix for any small e = — jp, the corresponding unique smooth g 6 ^ £ H°°(M. 3 ) stationary 



state to (1.11) with zero mean velocity and temperature fixed by the initial data. Then, we can 



show the following result: 

4.1 THEOREM. Given M the Maxwellian with zero mean velocity and temperature given by the 
initial temperature of fo, then there exist a positive constant C such that 



1 + I log e 



1/2 



for any e > small enough. 

Proof: Let g £ (t) be the solution to ( |1.11[ ) with initial data M, then 

||<& - M|| x i (R3) < \\g £ {t) - <4|| L i (K 3) + \\g £ {t) - M\\ L i 



Now, we are going to control each term separately. Since M G i^°°(M 3 ), then g £ (t) will satisfy due 
to Theorem II .31 that 

ll/(t)-^IUi(R3)<CE- r e-T / * 



for all t > 0. Here, we have made explicit the dependence on the restitution coefficient of the 
constants in the previous section. Actually, revising the discussion on the value of the constants in 
the previous section, one gets that 7' can be made as close to ^-7 as we want since our solution 



g £ (t) lies in ff°°(IR 3 ) due to Theorem 3.6 see last paragraph of the previous section. Moreover, 



we can fix e small enough and a close enough to 1 in such a way that 7 is as close as we want 
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to 2/15 in (3.1). For example, by choosing 1/11, then for a ~ 1 and for small enough e we have 
1/11 < j' < 16/165 = 2/15 • 8/11, where 1/11 is an arbitrary number. 

Concerning the behavior of the constants in front of the exponential in time function as e — > 0, 
it is not difficult, but tedious, to check that it does not degenerate to and is uniformly bounded as 
e — > in the case of the constants for the d% distance in Theorem 3.2 and in Theorem 3.3 However, 
the dependence on the restitution coefficient of the estimates in |13| Lemma 7.13, Proposition 7.30] 
leading to Lemma 3.4 degenerates as e —* as e~ r with an exponent r related to the regularity 
needed in the interpolation Proposition 3.8 One can estimate this degeneracy exponent r exactly 
depending on the regularity needed for having 1/11 < 7', but it is not important its exact value as 
we shall see below. 

On the other hand, we can use the Csiszar-Kullback inequality |14[ [20] together with the 
Logarithmic Sobolev inequality [191 122] to get: 

I(g%t))-I(M)> [ g ^ t ,v)ln^f^dv>hg%t) - Mf Ll 
J R s M [v) I 

Using Theorem |1.1[ we deduce 

U(t)-M\\ Lim < [2 (e^-l)/(M)] 1/2 

with 



2 + e + 15e 2 _ 1 + e _ 7 
4e 3 2~ ~ 2 

as e — ► 1. 

Finally, for suitable choice of a and for e small enough, we conclude 



I foe 



Ml 



L 1 



< 



2 e5 



et 



l)l(M) 



1/2 



+ Ce~ r e-ir 1 



for all t>0, and thus by Taylor's theorem, we get 



1 9lo 



Ml 



L 1 



< 



7ete2 £t I(M) 



1/2 



for all t>0. By choosing t = 11(1/2 + r)\ loge|, we obtain 

||^o -M\\ LHm < [77(1/2 + r)e I logele^^+^l 1082 ' J(M) ''"iCe 1 ' 2 
from which the announced result follows. 



1/2 



□ 



Appendix: The kinematics of inelastic collisions 

Here, we will review in detail the collision mechanism for inelastic collisions and the weak and 
strong formulation in two useful representations of the inelastic gain collision operator. We will 
perform in detail the relations between the collision frequencies in the different representations and 
for general interactions being of Maxwell type or not. Basically, these results make a summary of 
already known relations in particular cases written in [5] [18] but we believe this summary sets up 
the more general case once and for all. 
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A.l The kinematics of elastic collisions 



We begin by reviewing two ways of parameterizing the set of all elastic collisions in M 3 . The 
presentation has some unusual features that will be useful to our investigation of inelastic collisions. 
If particles with like masses and with velocities v and w collide elastically, so that both energy and 
momentum are conserved, then 

v + w 

and \v — w\ 



are both conserved. The conservation of the first quantity directly expresses the conservation of 
momentum (since the masses are the same), and then the conservation of the second follows from 
the conservation of energy and the parallelogram law: 



Therefore, let us introduce 



Note that 



If I 2 + \w\ 2 



v + w 
2 



v + w 



+ 



w 



v — w 



and 



u 
u 



z + 



-k and w = z — k . 



(A.l) 



(A.2) 



2 2 
Now consider an elastic collision (v,w) — > (y',w') where v' and w' are the post collisions 
velocities of the two particles. We define z', u' and k' in terms of the post collisions velocities v' 



and w' just as we defined z, u and k in terms of v and w in (A.l ). By the conservation of z and \u 



and (A.2), we have that 



and 



ii- 



iv 



k' 



(A.3) 



That is, by the conservation laws, only k changes, and the outcome of the collision is entirely 
specified by giving the change in the unit vector k — ► k', together with the initial velocities v and 



w. Hence the space of kinematically possible collisions is the set 



x S 2 with generic 



point (v,w,uj). The vector u is called the collision vector, and it is the additional parameter, 
beyond v and w, needed to specify the post collisional velocities v' and w'. This specification is 
then described by a bijective map C from 3 onto itself with C(v,w,u) = (v', w', u'), where v 1 and 
vJ are the post collisional velocities of the two particles, and oo and lo' are the collision vectors that 
are needed to determine k! from v and w, or in reverse, to determine k from v' and w' . 

The map C is called the collision map. There are two collision maps that are particularly useful 
for our purposes: the swapping map and the reflection map. The first has many mathematical 
advantages, due to its simplicity, but the latter has a closer connection with the physics of the 
collision process, and we shall need them both. 



A. 1.1 Reflection map: 

Consider a collision of two identical hard spheres in the center of momentum frame, so that z = 0. 
Let n be the unit vector pointing from the center of one particle to the center of the other at the 
moment of collision. (It does not matter from which to which; our expressions will be quadratic 



20 



in re.) The result of the collision is exactly as if each particle undergoes specular reflection upon 
striking the plane through the point of contact that has unit normal n. As a result, the velocities 
of both particles are reflected about this plane, as shown in Figure [TJ Clearly, the relative velocity 




Figure 1: Profile of the plane of reflection, with some normal lines indicated, along which the pre 
and post collisional velocities line up. 

u = v — w of the two particles is also reflected about this plane, so that we have v! = u — 2{u ■ n)n. 
This remains true if we translate out of the center of momentum frame. Also, since reflection is an 
isometry, k' = k — 2{k ■ n)n. We therefore let n stand in the role of uj, and we have 



k = k — 2{k ■ n)n 



and 



n 



n 



(A.4) 



From this and (A. 3) we obtain the full collision map: 



w 



n 



v — ((v — w) ■ n)n 
w + ((d — w) ■ n)n 



(A.5) 



n 



Let us use C r to denote the collision map C r (v,w,n) = (v',w',n') based on (A.5). Of course, n 
and — n lead to the same reflection, so that if we hold v and w fixed and vary n, we get a double 
cover of the set of possible post collisional velocities (v',w') resulting from (v,w). However, as is 
clear from the third line, as a transformation from H onto itself, C r is a bijection. 

The next thing to observe is that C r is a measure preserving transformation on S. Indeed, since 
the measure dvdwdn on H is also the measure dzdudn, going to polar coordinates in u, we have 
dv dw dn = dz d\u\ dk dn. Here dra and dk are both the uniform measures on S 2 with total mass 
4-7T, which is consistent with our use of polar coordinates. Since it is clear from Fubini's Theorem 
that dkdn = dk'dn', we have that dv'dw'dn' = dvdwdn. This gives us the collision kernel in the 
reflection representation 



(Q(f,f),<p) 



l 

4ir 



f(v)f(w)[(p(v') - <p(v)]&(\u\)B(k ■ n)dndvdw 



(A.6) 
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Here, 3> and B takes into account the difference rates of collisions and its dependence with respect 
to the strength of the relative velocity and the collision angle respectively and ip is any test function. 

Finally, since reflections are their own inverses, one sees that C r is its own inverse. Therefore, we 
may regard (v',w') as pre collisional velocities that may result in the pair (v, w) of post collisional 
velocities. Thus, in the "gain term" part of the integral for Q(f, /), one can change variables 



f(v)f(w)ip(v')B(k ■ n)dndvdw = 

Now, recall that B is even, and note that k' -v! = 
preserving property, we obtain 

f(v')f(w%(v)$(\u'\)B(k'-n')dn'dv'dw' 

is 2 

Thus, we have 



f(v')f(w')(p(v)&(\u\)B(k' ■ n')dn'dv'dw' . 

\u\. Using this and the measure 



k ■ n and I v! 





S 2 



f(v')f(w')(p(v)(f>(\u\)B(k-n)dndvdw 



(Q(f,f),<p) 



1 

47T 



S 2 



[f(v')f(w') - f(v)f(w)]cp(v)$(\u\)B(k ■ n)dndvdw 



which allows us to write the collision kernel in the strong form: 

1 



Q(f,f)(v) 



[f(v')f(w') - f(v)f(w)M\u\)B(k • n)dndw . 



(A.7) 



4vr j R 3 j S 2 

Note that to specify everything needed to compute Q(f,f), one needs only the first two lines of 



(A. 5), and this is all that is usually written down in discussions of elastic collisions. However, in 



the inelastic case, the collisions are not reversible, since they dissipate energy, and so there is not 
such a simple relation between the pre and post collisional velocities. Thus, more care is required 
in the passage from the weak form of the collision kernel to the strong form, and it will be helpful 



to keep all three lines in (A.5), and remember that C s is a measure preserving bijection of E onto 



itself. Indeed, note that for fixed n, the map (v, w) i— > (v ' , w') is a two to one, so dropping the third 
line, we would not have a bijection. 



A. 1.2 Swapping map: 

As we have mentioned, there is another collision map which leads to a different way of writing 
down the collision kernel. This other way is less directly connected with the physics of hard sphere 
collisions, but it does have considerable mathematical advantages due to its simplicity, which is 
based on a simple swapping of ui and k. In this context, it is traditional to write a in place of u, 
and the very simple rule for computing k' and a' in terms of v , w and a is simply 



k' 



a 



and 



A: . 



(A.? 



As with (A. 4), we use this and (A. 3) to obtain the corresponding collision map: 

. v + w \v — w\ 



2 

V + w 



+ 



2 

v — w\ 



-a 



a 



-a 



(A.9) 



v — w 



\v — w\ 
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Let us use C s to denote the swapping map C s (v,w,a) = (v',w',a'). Like C r , C s is a measure 
preserving transformation on 5. As before, we note that the measure dvdwda on H is also the 
measure dzduda, and going to polar coordinates in u, we have dv dw da = dz d\u\ dk da. In this 
form it is clear that the collision map for swapping is a measure preserving map, dv'dw'da' = 
dvdwda. This gives us the collision kernel in the swapping representation: 



(Q(f, /),</>) 



l 

47T 



S 2 



f(v)f(w)[ip(v') - (p(v)]®(\u\)B(k ■ a)dadvdw 



(A.10) 



where <p is any test function, and as above k denotes the unit vector in the direction of v — w, and 
<]? and B gives the relative rates of the various kinematically possible collisions with respect to the 
strength of the relative velocity and the collision angle, and da is the uniform measure on S 2 with 
total mass 4-7T. 

Note that like C r , C s is its own inverse, simply because the map in (A. 8) is its own inverse. 



Thus, we may once more regard (V, w') as a pair of pre collisional velocities. Changing variables, 
and using k' ■ a' = k ■ a, and the measure preserving property, we obtain the strong form 



Q(fJ)(v) 



l 

4-7T 



M 3 JS 2 



[f{v')f{w') - f{v)f{w)]B{k ■ a)dadw 



(A.ll) 



Note that only the first two lines of (A. 9) are required to compute Q(f, /), and this is all that 



is usually written down in discussions of elastic collisions. However, it is worth noting explicitly 
that if one holds a fixed, and then just considers the transformation (v,w) i— ► (v',w') described by 
the first two lines, this transformation is not onto and not injective: While the direction of u is 
arbitrary, the direction of v! is always that of a. In the inelastic case, everything will be clear if we 
always keep in mind that the collision map is a bijection from H onto itself, and keep all three lines 
needed to specify this map. 



There is one more important point to notice: The integral in (A. 10) is unchanged if we swap v 
and w; i.e., make the change of variables (v, w, a) *— > (w, v, a). This transformation does not affect 
v' but it reverses the sign on k, and so we may replace B{k ■ a) in (A. 10) with B{—k ■ a) without 



affecting (Q + (/ ', /) , tp) . That is, only the symmetric part of B contributes to the collision kernel, 
and we may freely require, by symmetrizing if need be, that B is a symmetric function on [—1,1]. 
We shall always impose this symmetry requirement on our rate functions B. 



A. 1.3 Relations between representations 

Our final business in this section is to relate these two representation of the gain term, and to 



determine in particular the relation between <£, B and <£, B in (A. 10) and (A. 6). The crucial fact 
is that both maps (y, w, n) i— > (V, w' , n') and (v, w, a) t— > («', w' , a') yield the same pair (V, w') if a 
and n are related through a = k — 2{k ■ n)n; i.e., if a is the reflection of k about the plane normal 
to n. Indeed, since, for example, 



v + w v 
+ — 



w 



((v — w) ■ n)n 



((d — w) ■ n)n 



the first lines of (A. 9) and (A. 5) coincide if we relate n and a through 

[u ■ n)n = - — -a or, equivalently k — 2{k ■ n)n 



a 



(A.12) 



(A.13) 
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Since a is the reflection of k about the plane orthogonal to n, we recover n (up to a sign) in terms 
of a and k through: 



u — \u\a 



n 



or, equivalently 



k — a 



n 



(A.14) 



\u — \u\a\ \k — a\ 

The fact that we only recover n up to a sign does not matter; it is not n itself, but only the plane 
normal to n that matters in all computations we shall make. Finally, from (A. 13) doing the k- 
operation over the formula for k, we have 



\k ■ n\ 



1-k-a 



(A.15) 



These formulas will be very useful in relating the n representation and the a representation for 
inelastic collisions. 

To work out the relation between $, B and $, B, we need just one more identity [3j S]: For 
any test functions ^ on I 3 , 



S 2 



u — \u\a 



da 



2\u ■ n\ 



s- 



if ((n • n)n) dn . 



(A.16) 



To see this, observe that if we define 6 = cos (k ■ a) and % = cos l {k -n), then from the reflection 
relation between k, n and a in ( |A.13 ), we see that for < x ^ 7r /2, 9 = n — 2x, so that 

sm(6)d8 = 2sin(2x)dx = 4 cos(x) sin(x)dx = 4\k ■ n|sin(x)dx • 



Thus, using (A.15), 



da = 4\k ■ n\dn 



and 



dn 



4V 1-k-a 



-da . 



(A.17) 



Finally, by taking into account that with k fixed, n i— > a is a two to one cover of S 2 , and since 
(u — \u\a)/2 = (u ■ n)n as noted in (A. 13), we obtain the identity ( |A.16 ). 

Now since v' can be expressed either in terms of (u — \u\)a/2 or (n • n)n, by equating the right 
hand sides of (A. 10) and (A. 6) which requires $ = $ and B{k ■ a)da = 2B{k ■ n)dn since n i— > a is 
a two to one cover of S 2 and B is symmetric, we obtain from (A.17) that 



B(k ■ a) 



B(k-n) 



2\k ■ n\ 

Then since from (A. 13) we have k ■ a = 1 — 2{k ■ n) 2 , we obtain 

B U(l - *)/2 



(A.18) 



B{s) 

2V(l-*)/2 

where — 1 < s < 1 and — 1 < t < 1. 



and B(t) = 2\t\B(l -2t 2 ) 



(A.19) 



24 



A. 2 The kinematics of inelastic collisions 

The kinematics of inelastic collisions is easiest to describe starting from the n representation. As 
in the elastic case, viewed in the center of momentum frame, the collision behaves like a collision 
with a wall running normal to the direction vector n pointing from the center of one ball to the 
center of the other ball at the point of contact. However, in the case of an inelastic collision with 
a wall, the component of the reflected velocity normal to the wall is reduced, while the component 
parallel to the wall is unchanged. Thus the rule for updating the relative velocity after the collision 
is this: Let P n denote the orthogonal projection onto the span of n. Then 

P n (v' - w') = -eP n (v - w) and P^{v' - w') = P^{v - w) . 

Thus, our rule for updating the relative velocity after the collision is a sort of reduced reflection 

v' — w = (v — w) — (1 + e)((v — w) ■ n)n , 

where the restitution coefficient e satisfies < e < 1. If e = 1, then this is a reflection about the 
plane normal to n. If e = 0, this simply cancels out the component of v — w along n, which corre- 
sponds to a perfectly inelastic collision. This gives us the inelastic reflection map with restitution 




Figure 2: Profile of the plane of reflection, with some normal lines indicated, along which the pre 
and post collisional velocities line up in the inelastic case. 

coefficient e, denoted C r , e : C r>e (v,w,n) = (v',w',n') where 

1 + e 







V 


= v — 


< w' 


= w + 


n' 


= n 



2 

1 + e 



((v — w) ■ n)n 
((v — w) ■ n)n 



(A.20) 



Note that with v! = v' — w' , we have 



,'|2 



u 



M 2 + (e 2 -l)(u-n) 2 . 



(A.21) 
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and 



k' 



k — (1 + e){k ■ n)n 



(A.22) 



Vl + (e 2 -l)(fc-n) 2 ' 

We now have what we need to write down the gain term in the collision kernel in weak form 
using the inelastic reflection map: 



{Qt(f,f)M 



1 

4-7T 



s 2 



f(v)f(w)ip(v')$>(\u\)B(k ■ n)dndvdw 



(A.23) 



where <p and k are as before, and again, $ and B gives the relative rates of the various kinematically 
possible collisions, and dn is the uniform measure on S 2 with total mass Air. 

All formulas involving the Fourier transform, of which we shall make extensive use, are simplest 
in the u-parameterization, and so we must translate the above inelastic parameterization into terms 
of a. This is easy to do using (A. 13) and (A. 14). For example, the first line of (A. 20) can be written 
as 



. v + w u 1 + e v + w 1 — 

v = 1 (u ■ n)n = 1 

2 2 2 y ' 2 4 



e 1 + e. 
—u H — \u\a 



where (A. 13) was used to eliminate n in favor of a. In the same way, we translate the second line, 
— — i±£|«|cr, and it follows that v! = + ^|u|cj, so that 



and find w' 



and 



k' 



lu n — + — k .a 

(1 - e)k + (1 + e)a 
^2(l + e 2 )+2(l-e 2 )A;- ( T ' 



(A.24) 
(A.25) 



Finally, from (A. 20) we have (n • u')n = — e(n • u)n. Then from (A. 13) we have 
u — \u\a = 2{u ■ n)n = — 2e(u • n)n = —e(u — \u\a) . 

solving for a' we find 



a' = k' + e— r(k — a) . 
\u'\ 

Using the expressions just derived for and k' , we find 

(1 + e)k + (1 - e)a 



a 



v / 2(l + e 2 )+2(l-e 2 )A;-cj 



(A.26) 



This gives us the inelastic swapping map with restitution coefficient e, denoted C S)G : 
C Ste (v,w,a) = (v',w',a') where 



v + w 1 — e . 

— ^— + -^-( v ~ w > 

v + w 1 — e 



1 + e, 



\v — w\a 



w 



a 



( ^ 1 + e l 
[v — w) H ; — \v — w\a 



2 4 v ' 4 

(1 + e)fe + (1 - e)a 
^2(1 + e 2 ) + 2{1 - e 2 )k ■ a 



(A.27) 
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We now have what we need to write down the gain term in the collision kernel in weak form 
for inelastic collision in the a representation: The basic expression that defines the gain term is 

(Qt,<p) = ^[ I I f(v)f(wMv'M\u\)B(k-a)dadvdw (A.28) 

47 T 7R3 J R 3 J S 2 

where cp is any test function, and as above k denotes the unit vector in the direction of v — w, 
and $ and B gives the relative rates of the various kinematically possible collisions, and da is the 
uniform measure on S 2 with total mass 4tt. Let us remark that being the n-representation more 
physically meaningful, the u-representation is a nice mathematical device to derive estimates in 
both the weak and the strong formulation. 



A. 3 Strong form of the inelastic collision kernel 



Finally, we want to derive the strong form of the gain term, and for this we need the precollisional 
velocities v* and w*. That is, given a collision map C, we define (v*,w*,lu*) = C -1 ^, w, uj). It is 
very easy to invert the the transformation {v, w, n) i— ► (V , w', n') in the reflection parameterization 
(A. 20): Simply use a restitution coefficient of 1/e. This gives us C~^(v,w,n) = (v* , w* , n* ) where 

1 + e. 



( 




V* 


= v — 


< w* 


= w + 


n* 


= n 



2e 

1 + e 



2e 



((v — w) ■ n)n 
((v — w) ■ n)n 



(A.29) 



Note that with u* = v* — w*, we have 



\u*\ 2 = \u\ 2 + 



1 



1 (u-nf 



u 



1 + 



1 (k-n)'- 



and u* ■ n = —^-u ■ n. Combining these we have 



k* ■ n 



k ■ n 



^e 2 + (1 - e 2 )(k ■ n) 2 ' 



(A.30) 



(A.31) 



Now, the determinant of the Jacobian matrix for the transformation in (A. 20) is easily seen to be 

2 



1 - e 



1 + e 



independent of n, so that the Jacobian itself is e. Thus, the determinant of the Jacobian matrix of 
the inverse transformation is —1/e, also independent of n, so the Jacobian itself is 1/e. We thus 
can rewrite (A. 23) as 



{oi, <p) 



4ir 



f(v*)f(w*)<p(v)<S>(\u*\) B ( k *' U ^ dndvdw 

S 2 e 



(A.32) 



where ip and B are as before, and k* = v* — w* . Taking into account (A.30) and (A.31) and the 
fact that B is even, 

1 



Qt(f,f)(v) 
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S 3 



f{v*)f{w*)$J{\u\,k ■ n)BJ (k ■ n) dndw 



(A.33) 
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where Bf 



and 



e V ' \^/e 2 + {l-e 2 )s 2 J e 
5+(r,s) = $ (J-^e 2 + (1 - e 2 )s 2 ) . 



(A.34) 



(A.35) 



Then, using the identities in (A. 13) once more, we may translate (A. 29) into the a parameteri 



zation, and obtain C s ^(v,w,a) = (v* ,u>* , a*) where 

v + w 1 - e 



w 



2 4e 

v + w 1 — e 
+ 



(v — w H ; — K> — w\(j 

4e 

r ^ 1 + e l l 
[v — w) ; — \v — w\a 



a 



2 4e v 4e 

(l + e)fc- (l-e)a 
" V2(l + e 2 )-2(l-e 2 )/fc-a 

Along the way, we find, using ( A.13| ) and ( A.14[ ) as before, 



(A.36) 



\u 



*|2 



U 



and 



Combining these we have 



|fc* - T»| 



2(2 ((l + e 2 )-(l-e 2 )(k.a)) 

(1 + e)a - (1 - e)k 
\/2(l + e 2 ) — 2(1 — e 2 )k ■ a ' 

' 1-k-a 
(1 + e 2 ) — (1 — e 2 )k ■ a ' 



(A.37) 
(A.38) 

(A.39) 



Now, we use (A. 17) and (A.37) to translate B{k* ■ n)dn into terms of a. We find, remembering 
that B is even: 



B{k * • n)dU = ¥ ( V (l + e 2 )-(l-e 2 )k.a V 



Next, we use (A. 19) to express this in terms of B in place of B: We obtain: 

(l + e 2 )k-a- (1 -e 2 Y 



B(k* ■ n)dn = B 



V2 



(1 + e 2 ) - (1 - e 2 )k ■ a J ^(l + e 2 ) - (1 - e 2 )k ■ a 



da . 



Moreover, using (A.31), we have 



*(|u*|) = *(|«*|) = $ (^L^{l + e 2 )-{l-e 2 ){k.a) 

Therefore, defining the function Bf by 

(l + e 2 )s-(l-e 2 ) 



B+(s) = B 



V2 



(1 + e 2 ) - (1 - e 2 )s J ^(l + e 2 )-(l-e 2 ) S e 
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and the function 



*+(r, a) = <J> (^V(l + e 2 )-(l-e 2 ) S ) 



we have 

Qt(fJ)(v) = - r / /(«W)^(l«U-^(fe'^- (A. 40) 

47r JlR3 J52 

The change of variable formulas that we have deduced in this section, suffice not only to allow 
us to write down the strong form of the collision kernel, but also to prove the following somewhat 
more general result: 

A.l THEOREM. Let K be any continuous real valued function on H x S. Then 

K[C~l(y, w, a), (v, w, o~)]$f(\u\, k ■ a)B^(k ■ a)dadwdv 

K[(v, w, a),C 8>e (v, w, a)]<f>{\u\)B(k ■ a)dadwdv (A.41) 



1 

47T 



s 2 

1 



4"7T 



s 2 



and 
1 

47T 



K[C r e (v, w, n), (v, w, n)]$g (\u\, k ■ n)B^{k ■ n)dndwdv 
s~ 

= 1-1 I I K[(v,w,n),C re (v,w,n)]$(\u\)B(k -n)dndwdv , (A.42) 
4vr J R 3 J R 3 J S 2 



where $ = $ and B and B are related by (A.19). 
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